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A More Improved Lattice Action for Heavy Quarks 
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We extend the Fermilab formalism for heavy quarks to develop a more improved action. We give results of 
matching calculations of the improvement couplings at tree level. Finally, we estimate the discretization errors 
associated with the new action. 



> 

o 
o 

m 
O 

> 

43 



X 



1. INTRODUCTION 

One source of uncertainty in lattice QCD comes 
from the discretization effects in numerical simu- 
lations. Here we extend the Fermilab formalism 
to higher dimensional operators to obtain a more 
improved action for heavy quarks, aiming to re- 
duce this uncertainty to 1-2%. 

The method starts with Symanzik £Q, who in- 
troduced a local effective Lagrangian 



Cos = £qcd + J2-j a 



dim C-j 



(1) 



of a continuum field theory. The Cj are higher- 
dimensional operators. The lattice action is ad- 
justed so that the short-distance Kj coefficients 
vanish, at least approximately. 

The Fermilab formalism 0] takes Wilson fer- 
mions with the clover action 3 and introduces 
different couplings for space-like and time-like in- 
teractions. Ref. 0] considers operators up to di- 
mension five. Ref. adds dimension-six oper- 
ators to this action. In this paper, we give a 
more complete action which has dimension six 
and some dimension seven operators and deter- 
mine the couplings by matching at tree level. 

2. QUARK BILINEARS 

In this section, we consider the quark bilinears. 
The dimension-five (-six) interactions are listed in 
Ref. 0] (Ref. Redundant directions are iden- 
tified by held transformations. They are used to 
eliminate the higher-order time derivatives, which 
would give large errors for large masses. We do 
not repeat this analysis here. 



The corrections in Eq. should be clas- 
sified by a power-counting scheme. For light 
quark physics the only relevant ratio is Aa. For 
heavy quark physics, however, the power count- 
ing should follow HQET for heavy-light systems 
or NRQCD for quarkonium. In HQET, operators 
are assigned a power of A = A/toq; in NRQCD, 
a power of relative velocity v. In these power- 
counting schemes, some dimension-seven interac- 
tions are commensurate with dimension-six inter- 
actions. As seen in the list in Table they scale 
as A 3 or v 6 . Further dimension-seven operators 
carry power A 4 (or higher) in HQET and v 8 (or 
higher) in NRQCD. 

Including the needed dimension-seven opera- 
tors, our lattice action takes the form 



S* 



So + Sb + Se 



^{~f-D,D 2 }Tjj + a 2 c 2 / ipjiDlip 



a 2 ci 



+ a 3 r 3 J i>{D 2 ) 2 ^ + a 3 r 4 J ipDji/j 
+ a 2 c 5 J ^{7 • B}ip 
+ a 2 c 6 / ip(D ■ E - E ■ D)ip 



+ a 2 c 7 J ipj ■ (D x B + B x D)ip 

+ a 2 c$ J ^{7 D,a E}i[> 

+ a 3 r 9 J ^{D 2 ,iV-B}4, 

+ a 3 r 10 [ ^^{i^B^D 2 }^ 
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Table 1 

Dimension six and seven bilinear interactions in 
C e g with HQET and NRQCD power counting. 



dim 


Operator 


HQET 


NRQCD 


6 




A 


V 4 




Ql ■ DQ 


A 


V 4 




Q{j D,a E}Q 


A^ 


4 




Q{f D,iE B}Q 


A' 5 


v° 




Q[L> 4 ,7 • E]Q 


A 4 


(i 

!T 




Q^4(D ■ E - E ■ D)Q 


A z 


A 

V 




Qj-(DxB + BxD)Q 


A 3 


V 6 


7 


Q{D 2 fQ 


A 3 


v 4 




QDfQ 


A 3 


V 4 




Q{D 2 ,iE-B}Q 


A 3 


V 6 




Qj ■ DiH ■ B~f ■ DQ 


A 3 


V 6 




Q{D 2 j ,iY, i B i }Q 


A 3 


V 6 




Ei^j QiXiDjBiDjQ 


A 3 


V 6 




QDii'E ■ BD t Q 


A 3 


V 6 




QD ■ (B x D)Q 


A 3 


V 6 




QE-EQ 


A 3 


V 6 




Qi'E ■ (E x E)Q 


A 3 


V 6 




QBBQ 


A 3 


V 8 




QiE ■ (B x B)Q 


A 3 


V 8 



a 3 r 12 / ^7 • DiY, ■ B-y ■ Dip, 



(2) 



where conventional difference operators and 
clover fields are understood. We have included 
some terms that we expect to be redundant, to 
retain some flexibility to, say, speed up computer 
simulations. 

For brevity we do not introduce terms, corre- 
sponding to the last six rows of Table ^ The 
operators QDii'E ■ BD t Q and QD • (£? x D)Q 
are redundant. Matching of the last four rows in 
Table Q] requires a calculation of Compton scat- 
tering, which is is still in progress. 



3. MATCHING 

We derive matching conditions for the improve- 
ment couplings in Eq. (J2| by matching on-shell 
quantities for small momentum but any mga. We 
do this by calculating the energy momentum re- 
lation and the scattering of a heavy quark off a 
background field. The energy momentum relation 
has been studied 0] and summarized |Hj before. 
It gives us conditions on c\, C2, r 3 and r 4 . Re- 
sults for the chromoelectric terms ce, cq and cs 
are given in Ref . , determined from the tempo- 
ral component of the background field. Scatter- 
ing from the spatial component of the background 
field yields results for the chromomagnetic terms. 
Our results yield the following matching condi- 
tions. The energy momentum relation gives: 

4C 4 (1 - C 2 ) 3r s 2 C 4 



16Cci 



[to (2 + to )] 2 
C 3 [2C + 4r s (l + m ) 



(1 +m ) 2 
■6r,C7(H 



mo)] 



to (2+to ) . 



mo (2 

vie 



too) 



3a3 



8r 3 



(3) 



2(l+m ) L (l + m ) (l + m ) 2 
where r 3 can be chosen to be zero. Also 

0^ 2 + 6 C c 2 + (, s C-24, 4 )^±^l 

8(1 + mo) 

We do not set r 4 = or r s <^/24 but choose it 
below so that rn = 0. For the chromoelectric 
interactions one finds 



(4) 



C 2 4 



c 8 - 



2 m (2 + m ) = C 2 (C 2 -1) 

(1 + TOo) 



mo(2+TO ) (1+mo) 
r 2 TO (2 + toq) 
4(1+ too) 2 ' 



c 6 = 0. 



(5) 
(6) 



For the chromomagnetic interactions one finds 
c B = r s , c 7 = 0, c 5 =ci, r 9 =r 3 , (7) 

!>'.., -■24/-,) - - —rj^. (8) 

(9) 
(10) 

In summary, we may set eg = Cj = eg = r 3 = 
rg = rn = r\i = 0, leaving only five non-zero 
couplings: c\, C2, r^, C5 and r\o- 



n.2 



1 

~48' 
48 v s 

0. 



ci, 
1 

12' 

24r 4 ) - \c B C, 
6 
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4. ERROR ESTIMATES 

We now estimate the discretization effects. It 
is advantageous to use a heavy-quark description 
of cutoff effects UJ, because then contributions 
from, say, Qa ■ EQ and Q{"f ■ D. a ■ E}Q are 
automatically combined. Thus, we express the 
uncertainties as 

Error, = (cj at - Cf ont ) (O,) (11) 

where Cj are HQET/NRQCD short-distance co- 
efficients for lattice gauge theory and continuum 
QCD Uj|. The mismatched coefficients can be ob- 
tained from the results of the previous section. 

The matrix elements (Oi) are estimated by us- 
ing HQET (NRQCD) power counting for heavy- 
light systems (quarkonia). Since higher-order 
calculations are not yet available we multiply 
the tree-level mismatch with a, to estimate the 
/—loop mismatch. This is a conservative estimate 
because the most pessimistic asymptotic behavior 
for the coefficients of the improvement operators 
is the same at higher orders as in the tree level 
formulas of previous section. To obtain the num- 
bers shown in Figure ^ we take a s = 0.25. 

Fig. n shows results for the chromomagnetic 
energy (labelled l/2ms), chromoelectric energy 
(labelled l/Am 2 E ), and the relativistic correction 
term (labelled l/8m|). See Refs. gE] for de- 
tails. The vertical lines represent a = 1/8 fm and 
a = 1/11 fm, corresponding to the unquenched 
MILC lattices |7]. The black (gray) curves rep- 
resent the charmed (bottom) quark. Note that, 
in most cases, bottom quarks have smaller dis- 
cretization effects than charmed quarks. Thus, 
Fig. G] shows concretely that the Fermilab formal- 
ism avoids errors of order (m,Qa) n . 

5. CONCLUSIONS 

We have presented a more improved action for 
heavy quarks with five new operators. Our error 
analysis shows that, to reduce errors below 1%, 
one-loop matching is essential for the coupling cb 
of the chromomagnetic clover term. Tree-level 
matching for the coupling of the chromoelectric 
term, ce, is usually enough, but one-loop may be 
necessary for charmonium. Tree-level matching 



HQET for heavy-light NRQCD for quarkonia 



relative error relative error 
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Figure 1. Truncation errors for the new lattice ac- 
tion. The dotted curves correspond to the unim- 
proved case. The dashed and solid curves show 
the errors from tree-level and one-loop matching, 
respectively. 

suffices for the other couplings in Eq. 
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